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ABSTRACT 


We extend previous studies of the physics of interstellar cloud collisions 
by beginning investigation of the role of magnetic helds through 2D 
magnetohydrodynamical (MHD) numerical simulations. In particular, we 
study head-on collisions between equal mass, mildly supersonic diffuse clouds 
similar to our previous study. Here we include a moderate magnetic held, 
corresponding to {3 = Vg/Vb = 4, and two limiting held geometries, with the held 
lines parallel (aligned) and perpendicular (transverse) to the colliding cloud 
motion. We explore both adiabatic and radiative {i] = Trad/Tcoii — 0.38) cases, 
and we simulate collisions between clouds evolved through prior motion in the 
intercloud medium. Then, in addition to the collision of evolved, identical clouds 
(symmetric cases), we also study collisions of initially identical clouds but with 
diherent evolutionary ages (asymmetric cases). 

Depending on their geometry, magnetic helds can signihcantly alter the 
outcome of the collisions compared to the hydrodynamic (HD) case. In the (i) 
aligned case, adiabatic collisions, like their HD counterparts, are very disruptive, 
independently of the cloud symmetry. However, when radiative processes are 
taken into account, partial coalescence takes place even in the asymmetric 
case, unlike the HD calculations. In the (ii) transverse case, the ehects of the 
magnetic held are even more dramatic, with remarkable diherences between 
unevolved and evolved clouds. Collisions between (initially adjacent) unevolved 
clouds are almost unahected by magnetic helds. However, the interaction with 
the magnetized intercloud gas during the pre-collision evolution produces a 
region of very high magnetic energy in front of the cloud. In collisions between 
evolved clouds with transverse held geometry, this region acts like a “bumper”, 
preventing direct contact between the clouds, and eventually reverses their 



motion. The “elasticity”, defined as the ratio of the final to the initial kinetic 
energy of each cloud, is about 0.5-0.6 in the cases we considered. This behavior 
is found both in adiabatic and radiative cases. 


Subject headings: ISM: clouds - kinematics and dynamics - magnetic fields ~ 
MHD; shock waves 
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1. Introduction 

Our understanding of the physical processes of the interstellar medium (ISM) of the 
Galaxy (and of external ones) has progressed tremendously observationally and theoretically 
in the last decade. Part of the required effort has been stimulated by possible implications 
for galaxy formation in the early epochs of the universe, but it is clear that many aspects of 
the subject pose specihc physics questions that are still unsolved and, therefore, interesting 
to study in their own right. 

That the ISM of our galaxy should present a multi-phase structure has been put 
forward for more than three decades, in its various versions including a two-, three- (and 
even four-) phase medium. The concept of a number of thermal phases coexisting in 
pressure equilibrium has now been developed further by Norman & Ferrara (1996) who 
found that, once turbulence is taken into account, a generalization to a continuum of phases 
is required. 

One of the aspects that has been recognized by essentially all the authors of the above 
mentioned studies as a crucial physical phenomenon in such a multi-phase environment 
is represented by “cloud collisions” (CCs). Of course, the term “cloud” might literally 
be appropriate only for an approach that is based on a somewhat simplihed thermal 
characterization of the ISM, whereas there is growing evidence that the dynamics of the 
gas, either in ordered or random/turbulent form, could govern its large-scale distribution of 
the gas. Nevertheless, collisions among fluid elements, in general, should be quite frequent 
and common in the ISM, and without sticking to any particular global model, the hrst 
aim of the present series of papers (Ricotti, Ferrara & Miniati 1997, RFM; Vietri, Ferrara 
& Miniati 1997; Miniati et al. 1997, Paper I) is to clarify the physics of such events, with 
special emphasis on the fate of “clouds” {i.e., , gas clumps) and the dissipation of their 
kinetic energy. 
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The main motivation for this study (and others) is to provide a hrm physical basis 
for global models of the ISM and galaxy evolution. Since the early paper by Oort (1954) 
several attempts (Field & Saslaw 1965, Habe et al. 1981, Struck-Marcell & Scalo 1984, 
Ikeuchi 1988, Vazquez & Scalo 1989, Theis et al. 1992, Jungwiert & Palous 1996) have been 
made to interpret the observed properties of galaxies (as, for example, their star formation 
history, ISM phase evolution and chemical evolution) as regulated by cloud self-interactions 
and with the environment (other phases, radiation held, gravitational potential). Basically, 
the assumption is that the ISM is mostly in cold gas clumps (he., clouds) with a spectrum 
of sizes, which are interacting, coalescing, forming stars (above a critical mass), fragmented 
again by energetic stellar events and re-injected into the ISM again. The model of Vazquez 
& Scalo (1989), for example, includes heuristic prescriptions on the fate of collisions through 
a function fc of the relative velocity. In this way they have been able to explain a number 
of interesting nonlinear behaviors as star formation bursts, whose properties are sensibly 
inhuenced by cloud interactions. Similar remarks can be done for models using cloud 
interaction physics to predict the phase interchange in the ISM (Habe et al. 1981, Ikeuchi 
1988). 

Several physical aspects of cloud collisions have been theoretically investigated in 
the past (Stone 1970a,b; Smith 1980; Hausman 1981; Gilden 1984; Lattanzio et al. 1985; 
Klein, McKee & Woods 1995); a detailed overview of the characteristics of CCs is given in 
Paper I and RFM. In spite of the fact that a Galactic collision event is difficult to observe 
because of its infrequency, short duration, emission and identideation, a growing amount of 
observational evidence for the events is slowly accumulating: some examples are found in 
NGG 1333 (Loren 1976), in Heiles cloud 2 (Little et al. 1978), in Draco (Rohlfs et al. 1989) 
and in NGP (Meyerdierks 1992). These detections often correspond to collisions involving 
clumps hosting a star formation region (Loren 1976). More recently Vallee (1995) has 
collected convincing data for a cloud collision event toward IRAS 2306-1-1451. 
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So far only minor attention has been devoted to magnetized CCs . Some pioneering 
analytical MHD work can be found in the literature (Clifford & Elmegreen 1983) but 
numerical studies have been overwhelmingly limited to hydrodynamical calculations. This 
is surprising, since by now magnetic helds have been detected throughout our galaxy by 
a number of dedicated experiments. Observations (Spitzer 1978, Zeldovich, Ruzmaikin, & 
Sokoloff 1983 and references therein) suggest that its orientation is mainly parallel to the 
galactic plane and, according to some authors, it becomes toroidal at high latitudes (Gomez 
De Castro, Pudritz, & Bastien 1997). The magnetic held is further believed to consist of 
a mean systematic component and of a random one. The strength of both is found to 
be approximately a few /xG. The evidence is provided by Faraday rotation, synchrotron 
radiation emitted by energetic electrons (cosmic rays), starlight polarization and Zeeman 
effect measurements (see e.g., Zeldovich, Ruzmaikin, & Sokoloff 1983, for more details). 
Direct information about the magnetic held along the line of sight (Ry) in galactic dihuse 
clouds, is obtained by observing the Zeeman splitting of the 21 cm radio line. With this 
technique the magnetic held strength is found to range on average between 3 and 12 pG 
both in H I and CO dihuse clouds (Myers et al. 1995). Also Heiles (1989) hnds Ry ~ 6.4 pG 
observing “morphologically distinct H I shells”. 


Recent measurements carried out by Myers & Khersonsky (1995) have substantially 
improved our knowledge of the properties of magnetic helds in interstellar clouds. For the 
H I dihuse clouds in their sample, \ogXe is found to range between -2.7 and -4.9, Xe being 
the electron fraction. The kinetic Reynolds number. Re {=vr/v), and even the magnetic 
Reynolds number, Reju {=vI/i'm) turn out to be large for n ~ a few km sec“^ and r ~ ~ 1 


pc. Those results validate our use of an ideal MHD code (see §^) for these simulations. 


Large Reynolds numbers are very familiar in astrophysics and characterize non-viscous 
hows. In addition when Re^ S> 1, the held lines are well coupled with the neutrals and the 
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magnetic flux is frozen into the fluid. Finally the ambipolar diffusion time is given by 

= .. yr. (1-1) 

and therefore the held should not decay through this process, in the timescales relevant 
for the clouds under study. In particular this ambipolar diffusion should affect the clouds 
neither during their propagation through the ISM prior to, nor during the collisions. 

Real structures in the ISM will have complex geometries, so any attempt to model 
specihc interactions in detail will require 3D simulations. Yet, ours is the hrst explicit MHD 
study of CC, where the physical effects due to the presence of a magnetic held are being 
investigated. Because any 3D studies will certainly require interpretation of very complex 
patterns and behaviors, we anticipate those works with an explorative 2D study that should 
contain many of the same physical behaviors, and, being far simpler to understand, ohers 
a practical basis for comparison. Several important and fundamentally diherent aspects of 
the physics of CCs are expected to come out of 3D calculations. One important example is 
the appearance of Kelvin-Helmholtz (KHI) and Rayleigh-Taylor (RTI) instabilities of the 
cloud surface along the cylinder axial direction, which are suppressed in 2D calculations. 
This might be particularly relevant when the magnetic held is transverse to the cloud 
motion, because in this case instabilities on the plane perpendicular to the cylinder axis 
can be suppressed by the development of an intense magnetic held at the cloud nose (Jones 
et al. 1996). This limits the time over which the 2D hows are really representative. At 
the same time, an initial look to our preliminary result of 3D single cloud calculations 
(Gregori et al. 1998) reveal that part of the MHD structure relevant to CCs developed by 
three-dimensional clouds is qualitatively similar to that seen in 2D clouds. Even though 
other quantitative diherences must occur (see §^), that result certainly supports the validity 
of our approach consisting of an initial explorative study of this yet uninvestigated problem. 


Since our objective is the examination of explicit MHD ehects in 2D cloud collisions. 



we try to follow as closely as possible the analogous HD simulations presented in Paper 
I. The plan of the paper is as follows. In Sec 2 we give some general considerations and 
introduce the relevant physical quantities of the problem; in Sec. 3 we briefly describe the 
numerical code and the experimental setup. Sec. 4 is devoted to the results, which are 
discussed and summarized in Sec 5. 


General Considerations and Parameters of the Problem 


2.1. Gas Dynamics 

In this section we review some basic aspects of purely hydrodynamical CC (HD CC), a 
problem already studied in great detail by previous authors {e.g., Paper I, Klein et al. 1995 
and references therein). The natural timescale for CC is given by 

Rr 


'^coll 


( 2 - 1 ) 


which is approximately the time required for the shock generated by the collision to 
propagate across the cloud radius, Rc- We suppose that clouds initially have a circular cross 
section. The main parameter determining the character of non self-gravitating HD CCs is 
(Klein et al. 1995) 


V = 




rad 


( 2 - 2 ) 


’^cRc 

where Wad = ^cVcTrad IS the radiative cooling column density through one cloud, is the 
cloud number density, and Trad is the cooling time (Spitzer 1978). Combining Eqs. |2-1| 
and |2-2| we have rj = Trad/Tcoii- Accordingly, if r/ < 1 signihcant radiative cooling takes 
place during the collision; when r; S> 1 emission processes become unimportant and the 
flow behaves adiabatically. In Paper I we concluded, after several low resolution tests, that 
the latter condition can be well represented by the weaker relation p > 1. Once the two 
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phase model for the ISM is assumed, for a given cloud velocity these conditions imply that 
impacts between larger clouds are more influenced by radiative cooling (RFM, Paper I). 

Since much of the basic physics of CCs has been explored using head-on events, which 
can provide a standard for comparison, in this paper we focus entirely on head-on CCs . In 
general, head-on symmetric HD CCs evolve through four main phases (Stone 1970a, 1970b); 
namely, compression, re-expansion, collapse and under some circumstances dispersal (Paper 
I). The occurrence of these four phases was first pointed out by Stone in his pioneering 
work (1970a, 1970b). It was also subsequently conhrmed by high resolution hydrodynamic 
calculations (Klein, McKee & Woods 1995, Paper I), which allows for the highest resolution. 
Similar results were also found through SPH calculations {e.g., Lattanzio et al. 1985, 
Lattanzio & Henriksen 1988). However, the limitations of these calculations, both in 
resolution and in the SPH method of simulation, have led to some misinterpretations of 
the physics of CC, like the hypothesis of “isothermality”, as discussed in Paper I. To a 
certain extent the same depiction of four phases can be drawn for magnetohydrodynamical 
CCs (MHD CCs) as well, although as we shall see, there are some important differences. 
However, we adopt this general terminology as a useful tool to refer to the various stages of 
the CCs . 


2.2. Cloud Propagation through a Magnetized ISM 

The HD of a dense cloud moving into a low density medium has been central to the 
stndy of several anthors (Jones et al. 1994, 1996; Schiano, Christiansen & Knerr 1995; 
Murray et al. 1995; Vietri et al. 1997; Malagoli, Bodo & Rosner 1996) who generally 
concentrated on the growth of KHI and RTI in snch conditions. We refer the interested 
reader to these works for an exhanstive description of this topic. When considering the 
motion of a clond through a magnetized medinm, new parameters, in addition to those 
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introduced so far (O. must be considered. In particular, the initial magnetic field 
is completely defined by its strength and orientation. The former parameter is usually 
expressed in terms of 



(2-3) 


where Pg and pb = are the gas and magnetic pressure, and M = Vc/cg and 

Ma = Vc/{B /A/dvrp) are the sonic and Alfvenic Mach number respectively. The magnetic 
held orientation is, in general, determined by two angles; in 2D simulations, with the held 
lying in the computational plane, they reduce to 9, the angle between the cloud velocity 
and the held lines. As long as the unperturbed magnetic held is dynamically unimportant 
(/3o 1) the initial evolution of a MHD cloud is similar to a HD one. This is the case we 

consider below, where we adopt (3 = 4. So, in the presence of a weak held, a stationary bow 
shock develops on a timescale Tbs ~ ‘^'Tcou- Further, a “crushing” shock is generated and 
propagates through the cloud with relative speed Vcs — ^ timescale {crushing 


time) 


2Rc 2RcX^ 


(2-4) 


where y = pdPi is the ratio of the cloud and intercloud medium densities. Finally a low 
pressure region (wake) forms at the rear of the cloud and, interacting with the converging 
how rehected oh the symmetry axis (X-axis), generates a relatively strong tail shock. 
However, over time, new features develop in response to the magnetic held. Jones et 
al. (1996), in a 2D study of individual MHD supersonic clouds, identihed several of these 
for an adiabatic, high Mach number cloud with modest density contrast, y = 10. Mac 
Low et al. (1994) also studied the MHD evolution of 2D individual, shocked clouds, which 
behave in a qualitatively similar manner. Utilizing these works, in the next paragraphs we 
give a brief review of the main features that result from inclusion of a magnetic held during 
the interaction with the intercloud medium. 


First, the orientation of the magnetic held with respect to the direction of the cloud 
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motion is particularly important. So far only two extreme cases have been published; 
namely, cloud motion parallel {aligned) or perpendicular (transverse) to the initially uniform 
held. We present elsewhere calculations with oblique magnetic helds (Miniati, Jones & 
Ryu, 1998). In the aligned case, the held lines, following the how, are swept over the cloud. 
As a result, those lines anchored at the cloud nose, are pulled, stretched and folded around 
the cloud. Eventually, these lines experience magnetic reconnection, forming new hux tubes 
passing around the cloud, somewhat like streamlines in a smooth how. In this region the 
magnetic held never becomes dynamically dominant, although its realignment around the 
cloud contributes to smoothing and, therefore, stabilizing the how. On the other hand, held 
lines are drawn into the cloud wake. Flow in the wake stretches lines anchored in the cloud 
material, causing the intensity of the magnetic held to increase. This feature, referred to as 
the post-cloud “hux rope” (Mac Low et al. 1994, Jones et al. 1996), is the only one that 
becomes magnetically dominated {(3 ^ 0.1) for an aligned geometry. A similar wake region 
forms also in the transverse held case, where the held lines drape over the cloud, converging 
in its wake. In this case, however, those held lines are anti-parallel across the symmetry axis. 
Above and below the symmetry axis the held structures in the wake initially are relatively 
uniform with a very sharp transition between them corresponding to a thin current sheet. 
Classically such a thin current sheet is unstable to the resistive, “tearing-mode” instability 
(e.g., Biskamp 1993, p. 73), in which the current sheet breaks into line currents and the 
magnetic held reconnects across the sheet. The instability condition is that the thickness of 
the sheet is much smaller than its width (e.g., Biskamp 1993, p. 152 for details). Indeed, we 
see these sheet transitions break up into a series of closed held loops that are the signature 
of this instability (e.g., Melrose 1986, p 151). This rapid modihcation of the magnetic held 
topology is often called “tearing-mode reconnection” (e.g., Melrose 1986) and typihes the 
reconnection that occurs in our simulations. A very clear illustration of the evolution of 
one such example is shown in Miniati et al. (1998, Figure 3). Because of this behavior 
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the magnetic field intensities in the transverse case wakes are lower than in the aligned 
case. On the other hand, the field lines in front of the cloud are compressed and, more 
importantly, stretched around the cloud nose. In this region, unlike the previous aligned 
field case, reconnection does not occur for these field lines. Therefore, the magnetic field 
becomes very intense (10“^ < /3 < 10“^) on the cloud nose, forming a magnetic shield. As 
pointed out by Jones et al. (1996), the main reason of magnetic energy enhancement is the 
stretching of the lines as these are swept up by the cloud. The timescale for the growth of 
the magnetic energy is given in eq. 9 of Jones et al. (1996), in the spirit of a first order 
approximation quantity, as r ~ VcRc^ = '^coii- Soon however, nonlinear effects become 
important and a more realistic timescale, as long as 2D approximation is valid, is provided 
by (Miniati, Jones & Ryu 1998) 

r ~ T,oii. (2-5) 

Further details on the development of the magnetic shield and on its dependence on cloud 
characteristics can be found in Miniati, Jones & Ryu (1998). It is important to notice that 
the magnetic shield acts to prevent the growth of KHIs and RTIs on the cloud surface. 

When radiative losses are included, as in the HD calculations, the thermal energy of 
the compressed gas is lost, reducing its pressure and allowing the cloud material to be 
compressed to very large densities. Particularly in the transverse field case, the cloud aspect 
ratio (length, x, to height, y) is highly increased by this effect. 


3. Numerical Setup 
3.1. The Code 

Our CCs simulations have been performed using an ideal MHD code, based on a 
second order accurate, conservative, explicit TVD method. Details of the MHD code are 
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described in Ryu & Jones (1995) and Ryu, Jones & Frank (1995); a brief description of 
some aspects concerning the inclusion of cooling and of a mass tracer is given in Appendix 
A. The V ■ B = 0 condition is maintained during the simulations by a scheme similar to the 
Constrained Transport (CT) scheme (Evans & Hawley 1988; Dai & Woodward 1997), which 
is reported in Ryu et al. (1998). We have used the 2D, Cartesian version of the code. The 
computational domain is on the xy plane, and the Z-components of velocity and magnetic 
held have been set to zero. 


3.2. Grid, Boundary Conditions and Tests 

In each CC simulation only the plane |/ > 0 is included in the computational box, and 
rehection symmetry is assumed across the X-axis. The length scale is chosen for each case so 
that Rc = 1.0 and the computational domain is adjusted to minimize boundary inhuences, 
as listed in Table |I[ Since the grid is Cartesian, our clouds are actually cylinders, with 
axes in the Z-direction. The top and right boundaries are always open. The left boundary 
is open in the asymmetric cases (3 and 4 of Table P and is rehective (as the bottom 
boundary) otherwise, when there is a mirror symmetry to the collision. With this latter 
choice we are allowed to use only half of the grid and reduce the computational time of the 
calculation. Only the highest resolution calculations, characterized by 50 zones across the 
initial cloud radius, are presented here. Lower resolution (25 zones per initial cloud radius) 
tests were also performed in order to check for consistent behavior. It turns out that beside 
inevitable quantitative differences, for each case we studied there is absolute consistency in 
the CC outcome. In particular it is apparent the qualitative agreement between low and 
high resolution calculations, in the density distribution and magnetic held structures that 
form out of the CC. 


- 14- 

3.3. Initial Conditions 

In this section we will discuss the initial conditions for our CCs. It is worth pointing 
out from the very beginning that, following the results of Paper I, most of the simulated 
CCs involved evolved clouds, i.e., clouds that have propagated through the ISM for about 
Ter before colliding. Initially, individual clouds have a circular cross section and uniform 
density and are in pressure equilibrium with a uniform background medium; the magnetic 
held is also assumed uniform throughout the domain. The relevant parameters, whose 
numerical values are given below, are the density contrast y, the Mach number M and the 
cloud radius Rc. The exact thermodynamic quantities characterizing the initial equilibrium 
state are not particularly important as their memory is lost soon after the beginning of 
the cloud evolution: for a supersonic motion, the thermal pressure of the shocked gas and 
the ram pressure of the how are dynamically far more important than the initial pressure 
balance. In addition, as already pointed out in j j2.2| , the cloud motion through the intercloud 
medium produces a variety of features which strongly alter the initial conhguration. As 
a result, despite the simplicity of the conditions at the onset of the cloud motion, before 
the collision takes place both the gas and the magnetic held are characterized by a rich 
structure. A comparison with the collision of two unevolved clouds is provided in the results 
section (EHI). to show the importance of considering prior cloud evolution in this study. 

The initial values of the above parameters are the same as in Paper I and are listed 
below. These values are inspired by the most recent observational studies and thought 
to be representative of the magnetized ISM. We assume a specihe heat ratio 7 = 5/3 
throughout. The inter-cloud medium has a density n, = 0.22 cm“^ and temperature 
Ti = 7400 K. Clouds are characterized by a density contrast y = nc/nj=100, so that 
the cloud density and temperature are ric = 22 cm“^ and = 74 K respectively and 

Trad ~ 3.7 X 10^ yr. The sound speed in the inter-cloud medium turns out Csi ~ 10 km s“^ 
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and the equilibrium thermal pressure for the ISM Peq/ks = 1628 K cm“^. When initially 
set in motion parallel to the X-axis, each cloud has a Mach number M = Vcjcsi = 1.5, and 
therefore Vc ~ 15 km s“^. Setting Rc = 0.4 pc, yields Tcou = Rc/vc ~ 2.6 x 10^ yr. Since this 
implyes r] = Trad/Tcoii ~ 1.4 > 1, in accord with these collisions behave adiabatically 
and therefore the cooling can be turned off. In the following collisions involving these 
clouds are referred to as adiabatic cases. On the other hand, setting Rc = 1.5 pc we 
have: Tcoii ~ 9.7 x 10^ yr and therefore rj ^ 0.38. These are the radiative cases. Table ^ 
summarizes these cloud characteristics. 

With a magnetic held included, three new parameters with respect to the HD case 
have to be determined: the held orientation and strength. The orientation is determined in 
2D by a single parameter, namely the angle between the initial cloud velocity and magnetic 
held. We explore two cases: magnetic held parallel (aligned case) and perpendicular 
(transverse case) to the cloud velocity. For the initial strength of the magnetic held, 
conveniently expressed by the parameter Pq (see Eq. |2-3|) , we assume Pq = 4, corresponding 
to B = 1.2 fiG. It could be argued that this value is somewhat smaller than what is usually 
observed. However, during the following cloud evolution, the held is stretched and amplihed 
and in several regions becomes energetically dominant (/3 -C 1 and B > l/iG, § 11.2|) . In the 
resulting conhguration, therefore, the magnetic held inhuence is not highly sensitive to this 
choice. 

Initially, the Jeans length of our typical dihuse cloud is \j ~ 29 pc3> Rc- Even though 
large density enhancements are produced during the compression phase in symmetric 
radiative collisions, Xj never becomes smaller than the vertical size of the clouds. For this 
reason we have neglected self-gravity throughout our calculations (see also Klein, McKee, 
& Woods 1995). Since we have concentrated on dihuse clouds, as opposed to molecular 
complexes, this approximation is justihed. 
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Similarly to Paper I, we consider collisions between both radiative and adiabatic 
clonds. However, for brevity, we discnss here only the most signihcant new resnlts for the 
two different cases. Table ^ snmmarizes the parameters of the collisions discnssed below. 
For Cases 1 and 2 in Table [|, the clond begins its independent evolntion at t 
where t = 0 corresponds to the instant when the bow shock of the two clonds hrst tonch. 
Analogously, for Cases 3 and 4 the two clouds are placed on the grid at t = 0 with their 
bow shock next to each other, after evolving for ^/ 2 'Tcr and respectively. In these four 
cases the initial magnetic held was aligned with the cloud motion, whereas for Cases 5-7 it 
was transverse. Case 5 is the only non-evolved calculation we present: by that we mean a 
uniform cloud of circular cross section placed on the grid in such a way that its boundary 
is only 2 zones from the (rehecting) Y-axis at t = 0.0. Finally in Cases 6 and 7 the cloud 
starts its evolution at t = —^/aTct and then is treated as in Cases 1 and 2 respectively. 
Animations of each simulation have been posted on our World Wide Web site at the 
University of Minnesota. 


4. Results 

4.1. Aligned Field (BX, Cases 1-4) 

In the aligned held case MHD CC show many similarities with HD CC. These are 
illustrated in Figure |l|, where the evolution of various MHD cloud integral properties 
are plotted as a function of time. The curves of kinetic and thermal energy as well as 
Ycm closely resemble those for HD clouds in Figures 8 and 9 of Paper I and the occurrence 
of all the phases characteristic of a HD CC (see j^ |2.1| and Paper I) is clearly seen. In general 
a stronger compression is generated when the clouds have mirror symmetry across the 
impact plane and only a weak reexpansion takes place in radiative cases, because most of 
the thermal energy is radiated away. The top right panel of Figure displays the total 
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magnetic energy. Since its variations are related to compression and/or stretching of the 
held lines, this quantity gives an approximate measure of the overall interaction between 
the magnetic held and the gas. In the adiabatic cases (solid and dotted lines) the large 
expansion undergone by the cloud gas produces both signihcant stretching and compression 
of the held lines. As a result, during the re-expansion phase the total magnetic energy 
increases of about 30%. On the other hand, in all radiative cases the total magnetic energy 
suhers only slight variations. 

We now begin specihc comparison, considering the adiabatic collision of two evolved 
identical clouds with aligned helds. 


4.1.1. Symmetric Cases 


The evolution of Case 1 is reported in Figs. |^a and ^b, representing the held line 
geometry (left panels) and density distribution (right panels) at four diherent times 
{t/Tcoii =4.5, 24, 48.75 and 75). As complementary quantitative information. Figs. ||, § and 
^ provide, for the same times, cuts along the primary axes of density, thermal and magnetic 
pressure, respectively. Lines there correspond to t = A.bTcoii (solid), t = 2 Atcoii (dotted), 
t = 48.75rco;; (dashed) and t = 75 tcou (dot-dashed). Cuts in the top panels are along the 
X-axis whereas those in the bottom ones are along the Y-axis. 


As already suggested above, the evolution of the CC is substantially unahected by 


the magnetic held. Indeed the CC goes through the four phases mentioned in §2.1| as in 
the HD case. For the compression and reexpansion phase there is a close quantitative 
correspondence, as it can be inferred from a comparison of Figs. ^ (solid and dot line) and 
Figure 2 of Paper I (although lines in the two hgures do not correspond exactly to the same 
time). This is not unexpected, though, since the compression and reexpansion phases are 
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dominated by the high pressure of the shocked gas. 

During the reexpansion phase a thin shell of dense (dp* < p < lOp*) cloud gas forms 
behind the reverse shock of the expanding material. A long hnger appears on the X-axis, 
due to the fact that the reexpansion hnds an easy way through the flux rope, where the 
density and the pressure are quite low (Mac Low et al. 1994, Jones et al. 1996). Also, in 
contrast to the HD case, the shell boundaries in Figure |^a (bottom left) are quite sharp. 
This difference has an important physical base. In fact, although the initial magnetic held 
has not sufficient strength to inhibit the onset of KHI (Chandrasekhar 1961), nevertheless 
it is able to reduce and eventually stop its growth. In fact, as the magnetic held lines, 
frozen in the gas, get stretched in the turbulent how, their strength is increased until during 
eventual reconnection they redesign the how pattern to a more stable conhguration. The 
criterion for this held dominance is that the local alfvenic Mach number falls to order unity 
or less (Chandrasekhar 1961, Frank et al. 1996, Jones et al. 1997). This is evinced by the 
presence of several held line loops on the external side of the gas shell (lower left panel in 
Figure |^a). Inside the shell the magnetic held intensity has severely dropped. According 
to Figure the magnetic energy density (dotted line) has been reduced with respect to 
its initial value (solid line) by a factor ranging from 10^ to 10^. This cannot be accounted 
for by expansion alone; indeed complex reconnection processes at the beginning of the 
reexpansion phase are responsible as well. 

The collapse phase (at t ~ IS.lbTcou) is chaotic and turbulent: in and around the cloud 
gas the density distribution is rather clumpy and the magnetic held has a tangled structure 
(top panels of Figure and dash lines in Figs. | and |). This situation will persist until 
the end of the simulation (dot-dash lines in Figs. ^ and It is during this phase that 
the magnetic held produces a qualitative change in the evolution of the CC. In Paper I we 
showed that during the HD collapse phase the reverse shock propagates toward the inner 
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region of the expanded cloud gas, whereas the external layer is shredded by KHIs and RTIs. 
At the end numerous hlaments £11 a region of about the same size and shape of the shell at 
its maximum extent. In the present case, on the other hand, the magnetic field lines which 
have been stretched by the vertical expansion of the gas shell, begin to relax, accelerating 
the gas at the top of the shell (near the Y-axis) toward the X-axis (top panels of Figure 
^). At the end of the collision {t = 75 tcoii) the initial cloud material is confined in a layer 
beneath the relaxed magnetic field lines, with a mean density p ~ 4pj (bottom panels of 
Figure |^) and is still laterally expanding. Within this new structure, with a thickness of 
several initial cloud radii, we find a clumpy density distribution (dot-dash lines in Figure 
and a weak, irregular magnetic field (dot-dash lines in Figure |^). Outside it, on the other 
hand, the magnetic field has the same initial configuration but greater strength, often twice 
as much as its initial value (dot-dash line in bottom panel of Figure 1). 

The analogous radiative case (Case 2) is illustrated in Figure ^ and its quantitative 
properties plotted in Figure |^. Again, an aligned magnetic field does not seem to influence 
the collision. As in the HD case efficient radiative cooling allows a much higher gas 
compression (solid line in top panel of Figure |^) and inhibits the strong reexpansion that 
would be driven by the high pressure of the shocked gas (solid line in bottom panel of 
Figure |^). During these phases the clouds coalesce, generating a well defined high density 
round structure about twice as large as the initial single cloud (this can be inferred from the 
position of the sudden drop in the dot line in top panel of Figure |^). No collapse phase ever 
happens. In addition a narrow jet of gas, characteristic also of radiative HD CC, is formed 
and extends along the Y-axis. As it propagates transverse to magnetic field lines a high 
density spot is created at the leading edge. However those features are usually unimportant 
because they only involve a negligible fraction of the total mass. Of more interest is instead 
the final fate of the mentioned structure. At the end of our simulation {t = ^Otcou) its edge 
near the X-axis is still expanding along the X-axis with a speed v ~ O.lCsj. According to 
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the three density cuts in the top panel of Figure the cloud edge (located at the sharp 
drop in each line) has been expanding at roughly the same speed (~ O.lCsj) throughout 
the evolution. Based on this velocity we can, therefore, estimate the time for the 
density of the new structure to drop by a factor a. Assuming that only one dimension of 
the cloud volume increases (at the speed of 0.1 Cgi) as a result of its expansion, we have 
Ta ~ aRc/{0.1Csi) = alOMTcoii = a 1.5 x 10® yr (we used Tcoh ~ 10® yr from §|3.0|). Since 
inside the new structure p ~ 40 — 80, a can be as large as 5-10, and Tq, comparable with the 
time between two cloud collisions. 


4.1.2. Asymmetric Cases 

As in Paper I we produce a simple asymmetric CC by colliding clouds with the same 
initial mass, velocity and radius, but evolved individually for a different time interval before 
they collide. By contrasting with mirror symmetric cases we can begin to see properties 
that are symmetry dependent. In the following cases (Case 3 and 4) the two clouds have 
been evolved for about ^/ 2 'Tcr and for Ter respectively (Table 1). 

The evolution of the adiabatic asymmetric CC resembles the analogous symmetric 
Case 1. Figure || shows the held lines (top panel) and the density distribution (bottom 
panel) at the end of the re-expansion phase {t = 22.5tcoii). We can identify a relatively 
dense shell (p ~ 5pi) with a large clump of gas on top of it as well as loops of magnetic held 
lines generated by reconnection events. Some new features appear, however, as a result 
of the broken symmetry. An example is the long tail of cloud gas on the right hand side 
of Figure due to the unbalanced momentum distributions in the two clouds along the 
X-axis. Nevertheless they are not so signiheant as to alter the general character of this CC 
with respect to Case 1. Therefore we expect that the collapse phase and the remaining 
following evolution of this case will not diher qualitatively from that of Case 1. We also 
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point out that as long as the asymmetry does not prevent the development of a reexpansion 
phase, Case 1 can be considered as qualitatively well representative of adiabatic MHD CC 
with magnetic held aligned to the initial cloud motion. 

The radiative case (Case 4), shown in Figs. |a, |b, presents new and interesting 
insights. The most crucial part of the evolution is represented by the cloud interaction 
during the compression phase. During that phase the older and more compact cloud 
(C2), moving from the right, attempts to plow through the other cloud (Cl) (Figure ii 
t = 6.75 tcoii)- Some of the features that are visible in Case 2, such as the vertical jet of gas, 
can also be recognize here if one carefully accounts for the distortions due to the asymmetry. 

However, we give particular attention to the new feature on the X-axis toward the 
left of Figure ^a (X~ —5Rc). This is a compact clump with density about 5 x lO^pi and 
velocity along the X-axis, Vx ^ — ^/ 2 Csi- At the end of the simulation {t = 22.5 tcoii) it has 
expanded and is still moving toward the left along the X-axis (Figure ^). At this time the 
“mass tracer” variables allow us to conclude that, despite the large prevalence of gas from 
the more compact cloud (C2), the new clump is the result of a partial coalescence of the 
two initial clouds. Its mass is about 10% larger than the initial mass of either cloud. Its 
density varies between ISp* and lOOpi and the velocity pattern suggests that it is expanding 
along the X-axis with v ~ Cgi/Q. By the same argument at the end of ^ |4.1.1| we can, 
therefore, conclude that for this case ~ a 9 x 10® yr. Again, before the cloud disperses 
in the background medium, the newly-formed clump is likely to undergo another CC. The 
remainder of Cl’s gas (feature to the right), moving transversely to the magnetic held, has 
formed a long hlamentary structure and has created a sharp cusp in the held lines (Figure 
|b). 


This result strongly dihers from the analogous purely HD calculation, where we found 
that the collision produced a large, low-density-contrast hlamentary structure eventually 
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fading into the background gas. 


4.2. Transverse Field (BY, Cases 5-8) 

We start this section by presenting the case of a CC between two unevolved clouds 
(Case 5). This calculation is mostly intended to provide a reference case when studying 
evolved CCs and, thus, to emphasize the importance of initial conditions in calculations of 
this type. 


4 -2.1. Unevolved Adiabatic Collision 

The results of this calculation (Case 5) are shown in Figure |T^, where two density 
images are superposed on held lines. The left panel captures the reexpansion phase at 
t = 7.5 tcoii] it closely resembles the analogous HD Case 1 of Paper I and no signihcant 
difference from Case 1 (of this paper), where the magnetic held was aligned to the motion, 
can be pointed out. This all means that no major role is played by the magnetic held, as it 
has been the case so far for all adiabatic CCs. The shell density has typical values (~ lOp*) 
and the thermal pressure drives the reexpansion. RTIs develop on the shell surface near 
the Y-axis. However, at the end of the collision a thick layer, qualitatively similar to that 
formed in Case 1 and 3 forms along the Y-axis. We also point out that the evolution of the 
total magnetic energy (dash line in bottom left panel) is reliable only up to t dTcoii, which 
marks the exit of the forward blast shock from the right boundary of the grid, along with 
consequent outflow of magnetic energy. The right panel image is from t = SOtcou- RTIs, 
aided by the formation of a slow, switch-off shock, which aligns the field with the expansion 
of the shell, have formed a large huger expanding almost parallel to the Y-axis. 
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4-2.2. Evolved Symmetric Cases 

As already mentioned in and discussed by Jones et al. (1996), an individual cloud 
moving transverse to a magnetized intercloud medium develops a region of strong magnetic 
field known as the magnetic shield. That feature dominates the collisions of such evolved 
clouds from the start of their encounter. Consequently, magnetic held effects dominate the 
evolution in Cases 6 and 7, in striking contrast to Case 5. Consequently, it is clear that 
meaningful simulations of CCs may depend on understanding the held geometry in their 
surroundings and on allowing self-consistent magnetic structures to evolve before collisions 
take place. 

Case 6 is shown in Figs. 0a and 0b, both including held line geometry and density 
images superposed on the velocity held. In addition, the solid lines in Figure show the 
usual time evolution of integral quantities. Peaks in the thermal energy and corresponding 
valleys in the 17 m curve, are signatures of the bow shock precompression and the “collision” 
undergone by the cloud . Note the simultaneous kinetic energy decrease and the magnetic 
energy increase, primarily due to cloud interaction with the magnetic held (Miniati et 
al. 1998) at the beginning of the simulation and to the “collision” event later on. 

As the cloud approaches the collision plane (Y-axis) in Figure p!2| a (bottom panel), the 
held lines are highly compressed (top panel) generating a strong repulsive force, opposite to 
the cloud motion. Eventually all the cloud kinetic energy is converted into magnetic form 
and stored as magnetic pressure. Then the cloud stops and, as the held lines reexpand, is 
accelerated backward and its motion reversed (Figure [T^). Therefore we can state that the 
magnetic shield acts almost like an elastic magnetic bumper. At the “apex” of this reversal 
phase both thermal and magnetic energy peak, whereas the kinetic energy obviously is 
about null (Figure ^). The latter, however, at the end of the simulation returns to about 
60% its initial value (Figure |^) and the cloud velocity is ~ Cgi- Also thermal and 
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magnetic energy are back to their initial values. Therefore the only effect of the collision is 
to dissipate part of the cloud kinetic energy with no other major consequences. 


The same qualitative result is also obtained in the radiative Case 7, shown in Figs. h3| a 
and |T^ and in Figure 0 (dot line). Whereas the kinetic and magnetic energy evolve as in 
Case 6, the thermal energy and Wm parameter have now similar qualitative behavior but 
much lower values (the scale for the dotted lines in the two bottom panels of Figure 0, is 
10 times smaller than for the solid ones); thus, the evolved cloud is mostly supported by 
magnetic pressure. Figs. |T^a and show the magnetic bumper effect for the radiative 
case as well. The hnal kinetic energy is about the same as in the previous case (50% of 
the initial value). However, the cloud gas has extensively spread out, especially in the 
X-direction, creating a dense elongated structure. Because of the high density reached by 
the cloud during the reversal phase, the resulting cloud shape is probably limited by the 
diffusivity of the code. Nevertheless, the hnal outcome is different from the analogous HD 
and previous aligned held cases and is more similar to Case 6. The cloud neither disperses 
nor coalesces; its structure is, however, strongly distorted and its kinetic energy partly 
conserved. 


5. Summary &; Discussion 

We have investigated the role of the magnetic held in CCs, through high resolution, 
fully MHD 2D numerical simulations. This paper is an extension of the gasdynamical study 
presented in Paper I. Our aim is to provide a hrst step toward the understanding of the 
physical role of magnetic helds in interstellar dihuse CCs . In particular, we have studied 
magnetic inhuences on: (i) the hnal fate of the cloud after the collision {i.e., dispersal, 
coagulation, shattering, hlamentation); (ii) the evolution of cloud kinetic energy and (iii) the 
ehects of CCs on the magnetic held structure in and around the clouds. These simulations 
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represent only an initial attempt to study a complex problem. To identify the most obvious 
and simplest behaviors we have restricted the geometrical freedom of the flow to 2D. For 
all behaviors, but particularly for (hi), our results need to be conhrmed by more thorough 
and extended 3D calculations. The main results can be summarized as follows: 

• Adiabatic, aligned held CCs are disruptive (as in the HD case), both for symmetric 
and asymmetric events. The remnant consists of an elongated structure of low 
magnetic energy in which cloud and background gas are mixed together. 

• Addition of an aligned held to radiative, symmetric CCs does not change the fact 
that they dissipate most of the cloud kinetic energy, which leads to almost complete 
coalescence of the two clouds. During asymmetric collisions of this type coagulation 
takes place as well, but the hnal structure has a mass only slightly (~ 10% for Case 
4) larger than either cloud initial mass; little alteration of the magnetic held line 
pattern is seen. This result is important, since purely HD asymmetric collisions of 
dihuse clouds have been shown to be highly disruptive (Paper I, Klein et al. 1995 and 
references therein). 

• In 2D motion of clouds moving transverse to the magnetic held leads to the formation 
of a magnetic shield in front of each cloud. When two evolved clouds of that kind run 
into each other, a magnetic shield may prevent direct collision from taking place. In 
our simulations the clouds remain separated by a magnetic barrier and bounce back 
with a fraction e of the initial kinetic energy. According to our results e ~ 0.5 — 0.6 
for both the adiabatic and radiative cases. This is probably an upper limit in more 
realistic situations, including 3D, and especially oh-axis collisions, since the magnetic 
bumper may be less developed and other degrees of freedom {e.g., rotation) are 
available.In addition, if a third dimension were included, after the collision of the 
clouds the compressed magnetic shield would partially reexpand perpendicular to the 
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direction of the initial motion, thus reducing e value further. 

Despite the caveats mentioned, much of the character represented in the last bullet 
may be independent of the symmetry of the collision. This was tested in part through a low 
resolution 2D numerical experiment of an asymmetric head-on, transverse-held collision. 
Results were consistent with those cited. In general, the magnetic shield is expected to work 
at some level for off-axis and largely asymmetric cases; transfer of momentum and angle 
scattering would occur in a similar way as in head-on collisions. When the magnetic held 
is aligned to the cloud motion, some of the arguments discussed in Paper I for off-axis HD 
CC should apply here as well. In particular off-axis CCs with an impact parameter h 
should be well represented by our asymmetric cases. Also for adiabatic cases, we expect 
that even for b ^ Rc, CCs should produce a reexpansion that is strong enough to disperse 
the clouds. The asymmetric radiative case results, however, caution from extending the HD 
results to radiative MHD CC when b is comparable to Rc- Those cases must be investigated 
in the future. Finally for b ~ 2Rc the collision should produce only minor perturbation to 
the clouds. 

A striking difference exists between the outcome of CCs when aligned and transverse 
held geometry are considered. Therefore, in order to model correctly CCs in the ISM it is 
important to understand the conditions for the formation of the magnetic bumper. Two 
points are crucially important in this regard: (i) the assumed cloud shape and (ii) the initial 
conhguration. 

As for the former, no 3D MHD individual supersonic cloud numerical simulation has 
been published so far. Especially near the nose of the clouds we should expect to see some 
kind if magnetic shield develop. On the other hand divergence of the how away from the 
nose and transverse to the prevailing held orientation should advect held lines away from 
the nose, thus limiting its development and extent. The importance of that ehect will 
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depend on the geometry of the cloud. A “pointed cloud” will have only a very limited 
shield. For example, Koide et al. (1996), studying the propagation of extragalactic jets 
through a medium with an oblique magnetic held, hud that the held lines distort to let 
the jet pass through. This ehect reduces the strength of the magnetic shield and could 
be important for spherical clouds. However, for cylindrical or hlamentary clouds, with 
axes transverse to the motion, held lines may be trapped at the front of the cloud long 
enough to play a major dynamical role. We point out that elongated clouds are not simply 
convenient to our 2D approximation. Rather, elongated shapes are expected for clouds that 
form in a magnetized environment, where the support provided by the magnetic pressure 
is anisotropic {e.g., Spitzer 1978). As already pointed out by Jones et al. (1996), realistic 
clouds show strong shape irregularities where the held lines can penetrate, be captured 
and, therefore, stretched to form some sort of magnetic bumper. But even if not so, when 
a cylindrical cloud moves in a transverse magnetic held, the motion of the gas along the 
cloud major axis, away from the stagnation point, is certainly slower at the center of the 
cloud where the stagnation point is located, than near the sides. Since the held lines are 
frozen into the gas, the central region of the cylindrical cloud is where the held lines are 
held longest. Therefore, an uneven magnetic tension is applied on the cloud and a bending 
of the latter is produced in the central region. As a result it becomes more difficult for 
the held lines to slip by the cloud enhancing the deformation of the cloud and increasing 
trapping of the held lines. 

Finally we have just begun a set of preliminary, low resolution numerical calculations 
to be presented in a subsequent paper (Gregori et al. 1998). From those calculations we 
can anticipate that a magnetic shield always forms and in a fashion qualitatively similar to 
that seen in 2D cases. Moreover, as expected, initially elongated clouds develop a stronger 
magnetic shield than “spherical” clouds. However, in the latter case the tension of the 
magnetic held lines wrapped around the cloud produces a strong deformation of the cloud 



shape, which grows strongly elongated transverse to the plane containing the held and the 
motion thus facilitating the formation of the magnetic shield. 

On the other hand the effect of the initial conhguration on individual 2D cloud 
evolution has been investigated by Miniati et al. (1998) who study the inhuence of the 
initial held orientation with respect to the cloud motion {6), of the cloud density contrast 
(x) and velocity (M) on the magnetic bumper formation. They hnd that, as long as 2D 
approximation is valid and 6 ^ 30°, the timescale for the formation of the magnetic bumper 
is of the order of r ~ Tcoh- Since clouds are slowed down by the ram pressure of 

the impinging hux on a timescale Tde ~ X'^coih they also concluded that magnetic bumpers 
are more likely to develop around high density contrast, low Mach number clouds. 

Another important 3D issue is the reexpansion of cloud gas in the direction 
perpendicular to the computational plane (along the Z-axis). This ehect is important 
because it could in principle modify our previous conclusions for the non disruptive cases. 
However, as it turns out this only sets a limit on the length of the cloud major axis 
for the adiabatic Case 6. In radiative cases lateral reexpansion involves only a small 
fraction of cloud mass, independent of the Y or Z direction (see § |4.1.1|) . On the other 
hand, in the adiabatic case, the rarefaction wave generating such reexpansion propagates 
at the sound speed of the postshock gas, whose temperature has been enhanced by a 
factor ~ M^. Therefore, Vexp ~ Mcgi- If the length of the cylinder is 0 = fhe 
reexpansion occurs on a timescale Texp — ^zl'Vexp ~ fJ^Rc/(Mcgi) = fiTcoii- Consequently, 
since Trad/Texp = Trad/l^'Tcoii = h/hi then if /i 3> 1 (£2 3> Rc for Cylindrical clouds), the cloud 
behaves as it did in 2D radiative cases discussed before. 

The aforementioned reasons justify 2D simulations as a valuable starting point for 
the more complex 3D MHD CCs. As already mentioned along the paper, however, in the 
latter we expect to observe new behaviors not seen in 2D calculations. We expect those to 
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be mostly related to differences in the evolntion of the magnetic held. Since the latter is 
dynamically dominant in transverse held cases, new behaviors will be more apparent there. 
For example, even when the magnetic shield forms in 3D as well, if its strength is much less 
than in 2D, then we may expect diherent results from those reported in this paper (more 
precisely in § [4.2|) . Nevertheless, instabilities and, in general, hows along the third direction, 
as well as quantitative diherences in cloud features developed in 2D and 3D, will certainly 
ahect the dynamics of the collision. We have begun a series of 3D MHD cloud simulations 
to address these complex issues more fully. 
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for stimulating discussions. The work by FM and TWJ is partially supported by NSF 
grants AST-9318959 and INT-9511654 and by the University of Minnesota Supercomputer 
Institute. The work by DR was supported in part by KOSEF through grants 975-0200-006-2 
and 981-0203-011-2. AF acknowledges hospitality of University of Minnesota where this 
work was started. 


Appendix A 

We discuss here some details of the treatment of the radiative losses and of the mass 
tracing routine. Radiative losses have been taken into account using the same approach 
as in Paper I, to which we refer for a detailed description. The radiative correction that 
we have applied is quasi-second order accurate (see, e.g., LeVeque 1997 for a general 
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discussion). Let us represent symbolically the “numerical” equation as 

Vtq = Vs{q), (5-1) 

where and Vg are the temporal operator and the spatial plus source operator respectively 
and q is the set of variables describing our system. A second order accurate splitting from 
time step n to n -|- 2 for the operator Vg would be (Strang 1968): 

SiL^LySi SiLyL^Si (5-2) 

where and Ly are the differential operator with respect to the A-coordinate and 
F-coordinate respectively and S is the operator representing general source terms. The 
subscript | means that the operator is applied for only half of a time step. Instead, we have 
used 

SiL^Ly S LyL^Si, (5-3) 

which assumes that Si Si = S. Hence the “quasi-second” order description. In addition 

2 2 

we have suppressed cooling inside the shock thickness. Indeed, in the physical shock layer 
the flow should be non-radiative, because the crossing time of the real shock thickness 
(artihcially spread out by the code) is much shorter than the cooling time. Also, since 
density and pressure are not accurate inside the shock, but only adjacent to it, radiative 
cooling could become artihcially large, reducing the performance of the code. This turns 
out of particular importance for the MHD calculation. The radiative cooling function we 
have used is identical to that in Paper I and is fully described in Ferrara & Field (1994). 

It includes cooling due to free-free emission, recombination lines, and collisional excitation 
lines as well as heating terms provided by collisional ionization and ionization by cosmic 
rays. For the cosmic ray ionization rate we adopt the value 2 x 10“^'^ s“^ as determined 
from observations by van Dishoeck & Black (1986). This and other rates can be found in 
Ferrara & Field (1994) and references therein. We neglect dust, and particularly PAHs, 
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photoelectric heating. This is certainly a rough approximation as long as an accurate model 
of the multiphase ISM is concerned {e.g., Wolhre et al. 1995). However, our aim here is 
to build a simple, albeit reasonable, model for the two-phase ISM and concentrate on the 
properties of collisions that do not depend drastically on the details of the multiphase 
structure. When pressure equilibrium is imposed, a two-phase (cloud -|- intercloud) ISM 
structure results. 

Finally a routine, based on van Leer’s second-order advection scheme (van Leer 1976), 
has been included to track the fraction of cloud material inside each grid cell. This quantity, 
referred to as “mass tracer” or “mass fraction” (Xu & Stone, 1995), is initially set to unity 
inside the cloud and zero elsewhere. The mass fraction allows us to discriminate between 
the different components in our simulations, which are the two clouds and the intercloud 
medium. In this way we can calculate various quantities of interest such as each cloud’s 
kinetic and thermal energy as well as Wm, the F-coordinate of the center of mass of each 
cloud (Jones et al. 1996). These are used in the analysis of our results. 
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Table 1. Summary of 2D-MHD Cloud Collisions Simulations 


CASE 

rjo 

Clouds Ages 

Cl C2 


B 

Grid-Size‘s 

End Time® 

{j~coll) 

1 

adiabatic 

0.75 

0.75 

3 

Bx 

15 X 15 

75.0 

2 

0.38 

0.75 

0.75 

3 

Bx 

15 X 10 

30.0 

3 

adiabatic 

0.5 

1.0 

3 

Bx 

30 X 15 

22.5 

4 

0.38 

0.5 

1.0 

3 

Bx 

30 X 10 

22.5 

5 

adiabatic 

0.0 

0.0 

3 

By 

10 X 20 

30.0 

6 

adiabatic 

0.75 

0.75 

3 

By 

20 X 7.5 

30.0 

7 

0.38 

0.75 

0.75 

3 

By 

20 X 7.5 

30.0 


^-All models use /3 = 4, 7 = 5/3, y = pdPi = 100, equilibrium pressure Peq/ks = 1628 K cm“^. Also, 
at equilibrium, we have Ti = 7400 K and rii = 0.22 cm“^ for the intercloud medium and Tc = 74 K and 
ric = 22 cm“^ inside the clouds. 

P — '^radj '^coll 

'^This is the relative Mach number for the cloud pair when the clouds are first set into motion. It is referred 
to the intercloud sound speed, Csi- 

'^The grid size is expressed in units of cloud radius. One cloud radius i?c=50 zones. 

®The end time is expressed in terms of collision time TcoU, and represents the total time from the beginning 


of the collision. 







Table 2. Clouds Characteristics'^ 



Rc 

'^coll - ^cf 

Ter = 2Rc^tVc 

T ^ 

Trad — 2 nA 


(pc) 

(yr) 

(yr) 

(yr) 

Adiabatic 

0.4 

2 .6x10^ 

b.OxlO'^ 

oo 

0.38 

1.5 

9.7x10^ 

2 .0x10® 

3.7x10^ 


^■866 note a in Table |[ 

’^See note b in Table ^ 

■^The actual value of the cooling time for the smaller cloud is also Trad = 3.7 x 10^ yr. However, as 
explained in the text, since r; > 1 the cooling does not affect the collision of these clouds and therefore it has 
been turned off during the simulations. For this reason we have set Trad = oo in the table. 






FIGURE CAPTIONS 


Fig. 1.— Panels showing the evolntion of kinetic energy (top left), thermal energy (bottom 
left) and Um coordinates (bottom right), of the colliding clonds and of the total magnetic 
energy within the grid (top right) as a fnnction of time, for the aligned held cases. For Cases 
1 and 2 the two colliding clonds are identical to each other, so only one is displayed. For the 
hrst three panels, solid lines correspond to Case 1, long dash to Case 2. Dot and short dash 
lines respectively to cloud Cl (left) and C2 (right) in Case 3. Dot short-dash and dot-long 
dash lines respectively to cloud Cl (left) and C2 (right) in Case 4. For the last panel (top 
right) solid lines refer to Case 1, dash line to Case 2, dot line to Case 3 and dot-long dash line 
to Case 4. The kinetic, thermal total magnetic energy in these plots have been normalized to 
the sum of the cloud energy (kinetic -|- thermal -|- magnetic) and the background magnetic 
energy at the beginning of the simulation. Therefore, with respect to Paper I the vertical 
scale of the kinetic and thermal energy of the cloud is reduced by a factor 1.18. 

Fig. 2.— From top to bottom, the four pairs of panels show field line geometry (left) and 
density distribution at t/Tcoii = 4.5, 24, 48.75 and 75. Coordinates are expressed in units of 
cloud radii. Re- Field lines are contours of the magnetic flux, and correspond to change in 
the latter by a factor 7. Density hgures are inverted grayscale images of /(p) = Top 
panels show the compression phase (t = 4.5rcozz) with held lines swept up by the vertical 
outhow at the symmetry plane of the collision. The following panel shows the re-expansion 
phase (t = 2 Atcou)i with the formation of long huger at the back of the cloud, a large 
spot near the symmetry axis and a thin shell-like structure. Numerous closed held loops 
witness the occurrence of reconnection processes. The next two panels capture the collapse 
{t = 48.75rcozz) of the large spot pushed by relaxing held lines. The last two panels correspond 
to the end of the collision {t = 75tcou) with the formation of a thick layer of cloud gas and 
low magnetic energy. 
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Fig. 3. — Log density cuts through the grid along the X-axis (y=0.3i?c, top panel) and along 
the Y-axis (x=0.1 bottom panel) for Case 1. Solid lines refer to t = A.bTcoiu dot lines to 
t = 24:Tcoii, dash lines to t = 4:8.75 tcoii and dot-dash lines to t = 75 tcoii- 

Fig. 4.— Same as Figure |^, but for the logarithm of thermal pressure. 

Fig. 5.— Same as Figure |^, but for the logarithm of magnetic energy. 

Fig. 6.— Field line geometry (left panel) and density distribution (right panel) for Case 2 
at t = SOtcou- Line contours and density images as in Figure 

Fig. 7.— Cuts along the Y-coordinate {y = O.lRc) for the log of density (top panel) and log 
of thermal pressure (bottom panel) for Case 2. Solid line corresponds tot = i.bTcoii, dot line 
to t = 17.2 tcoii and dash line to t = SOtcou- 

Fig. 8.— Field line geometry (top panel) and density distribution (bottom panel) for Case 
3 at t = 22.5tcoii- Line contours and density images as in Figure 

Fig. 9.— Field line geometry and density distribution for Case 4 at t = 6.75tcoii (a) and 
t = 22.7)Tcoii (b) respectively. Line contours and density images are created as in Figure 

Fig. 10.— Density distribution with superimposed held line geometry for Case 5 at 
t = 7.5 tcoii (left panel) and t = SOtcou (right panel). Line contours and density images 
are created as in Figure |^. 

Fig. 11.— Plots as in Figure |I]but for transverse held cases. Dash lines refer to Case 5, solid 
and dot lines to Cases 6 and 7 respectively. The latter two lines have been multiplied by a 
factor 10 and 100 respectively, in order to make them readable on the same vertical scale as 
the dashed line. 
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Fig. 12.— Field line geometry and density distribution for Case 6 at t = 12rcozz (a) and 
t = SOtcou (b) respectively. Note that, since the magnetic held is quite stronger than in 
previous cases, contours of the magnetic hux (for the held lines) correspond to change in the 
latter by a factor 14 (instead of 7 as before). Density images, on the other hand, are as in 
Figure ||. 


Fig. 13.— 
t = 22.5 tcoii 


Field line geometry and density distribution for Case 7 at f = 9.75tcoii 


(b) respectively. Line contours and density images as in Figure 


(a) and 
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